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Optical entropy bounds for metal nanoparticles immersed in nonlinear optical media and for 
nonlinear dielectric microdroplets on metal surfaces are calculated near the frequency of the surface 
plasmon resonance. Similar to the Bekenstein-Hawking result for the black hole entropy, the entropy 
bounds in nonlinear quantum nanooptics may be expressed as the ratios of the droplet perimeter 
(nanoparticle area) to the effective Planck length (effective Planck length squared). 



Linear and nonlinear nano-optics of surface plasmon 
polaritons (SPP)i£ attracts great deal of current atten- 
tion due to wide range of its applications in such diverse 
areas as optical nanolithography 3 , sub-diffraction-limitcd 
microscopy^, quantum computing 5 , optoelectronics de- 
vice integration on sub- micrometer scales^*!, biosensing 8 , 
etc. Many of these applications are possible due to 
the fact that the SPP wavelength in a given frequency 
range may be much shorter than the wavelength of free 
space photons. Such tight spatial localization of the SPP 
field gives rise to considerable field enhancement and 
hence, to many interesting nonlinear optical phenomena, 
which occur at very low light level aV . As the num- 
ber of SPP quanta involved in these effects is reduced, 
quantum mechanical effects play more and more impor- 
tant role in nonlinear nanooptics. Quantum sources of 
SPP or SPASERS 11 - have been theoretically proposed re- 
cently, while experimentalists move close to realization 
of such quantum nonlinear nano-optical devices^. Thus, 
reaching good understanding of nonlinear quantum nano- 
optics is a very important current goal. However, as we 
have pointed out recently^, this is a highly nontrivial 
theoretical task. Its complexity may be illustrated by re- 
calling the analogy between the nonlinear quantum nano- 
optics and the quantum gravity, which has been indicated 
mi 3 .. 

It is well known (see for example ref^ and the ref- 
erences therein) that the refractive index n = (e/i) 1 / 2 
plays the role of an effective metrics in optics. The 
Maxwell equations in a general curved space-time back- 
ground gij(x,t) are equivalent to the phenomenological 
Maxwell equations in the presence of a matter back- 
ground with the following spatial distribution of the di- 
electric and magnetic permeability tensors £ij(x, t) and 
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If we extend this analogy to nonlinear optics phenom- 
ena, we must conclude that the expression for the dielec- 
tric constant of a Kerr medium in nonlinear optics: 

e d = e« +4n X ^E 2 (2) 



(where E is the optical electric field, e' 1 ^ and are 
the linear dielectric constant and the third order non- 
linear susceptibility of the medium, and for the sake of 
simplicity we omit tensor indices of e^ 1 ' and X ) pl a Y s 
the same role as the role of the Einstein equation in gen- 
eral relativity: 

Rik = — J-Tik, (3) 

cr 

where k is the gravitational constant, and we assume 
T = 0. Both equations define the influence of the energy- 
momentum tensor on the effective metric. For evident 
reasons nonlinear optics cannot emulate gravity exactly. 
However, from comparison of equations (2) and (3) it 
is clear that x^ m nonlinear optics plays roughly the 
same role as the gravitational constant in eq.(3): positive 
X*- 3 ^ results in self-focusing or, in other words, attractive 
interaction of optical quanta. 

The spatio-temporal scale at which quantum metric 
fluctuations play an important role in gravitational ef- 
fects is characterized by the Planck length L p (see ref~ 
for a recent review on quantum gravity): 

L p = (^) 1/2 w 1.6 x 1(T 33 cto (4) 

& 

Due to extreme smallness of this parameter, effects of 
quantum gravity are out of experimental reach. On the 
other hand, the effective nonlinear-optical Planck scale 
originally introduced ini 3 - does not appear that far out 
of reach. In order to avoid undesirable wave-vector de- 
pendence in the original expression from ref^, let us 
introduce the effective Planck scale as 

L%" = (ftx (3) c) V4 (5) 

This expression is based on simple dimensional analy- 
sis. Assuming \^ ~ 10~ 10 esu, as in experiments with 
4BCMU polidyacetylene films on gold in refiS, we obtain 
L e J5 k, 2.3 nm. While this value is still too small for 
usual nonlinear optical experiments, it plays an impor- 
tant role in quantum nonlinear optics of surface plasmon 
polaritonsi 3 -. In fact, the effective two-dimensional (2D) 
Planck scale as seen by the SPPs having frequencies near 
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the frequency of the surface plasmon resonance 2 is con- 
siderably larger. In order to understand this, we must 
recall the basic facts about the dispersion law of SPP 
(more detailed discussion can be found u&i^). 

Let us consider SPPs which propagate along an inter- 
face between a metal film and a third-order nonlinear 
dielectric, which is located on both sides of the film. In 
such a case the dispersion law can be written asi 3 - 
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where e m (w) is the frequency-dependent dielectric con- 
stant of the metal. If we assume that e m = 1 — w p /uj 2 
is real (where oj p is the plasma frequency of the lossless 
metal) and that d is large, the dispersion law can be sim- 
plified as 
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This dispersion law approaches asymptotically uj sp = 
LOp/ (l + e^) 1 / 2 at large wave vectors. The latter frequency 
corresponds to the surface plasmon resonance. The dis- 
persion law (6) also diverges near oj sp in the general case 
of a lossy metal film if the metal film thickness d is small, 
so that the imaginary part of the term ±2e d e~ kd com- 
pensates the imaginary part of e m . In addition, pumping 
energy into the dielectric by optical or other means may 
create a dielectric with gain in a given frequency ranged, 
so that the imaginary part of e d may compensate losses 
in metal and k again diverges near ui sp . Very-short wave- 
length (in the 10-50 nm range) plasmons with frequencies 
near lo sp were observed in a number of near-field optical 
experiments^*^. 

It is clear from eqs.(6,7) that around u> sp (near e m (uj) = 
—e d ) both phase velocity c p i and group velocity of surface 
plasmons tend to zero, and the effective two-dimensional 
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diverges as seen by the plasmons: 
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In a similar fashion, the effective (as seen by SPPs) 

(3) 

two-dimensional nonlinear susceptibility "xXd °f the ^i- 
electric diverges near tv sp : small changes of the three- 
dimensional dielectric constant e d due to the Attx^E 2 
term in equation (2) are perceived as very large changes 
of the effective two-dimensional dielectric constant (8) by 
surface plasmon-polaritons: 



L$ = (hx&c) 1 '* (10) 
is much larger then L e J* from eq.(5). As evident from 

eq.(9), the numerical value oilfJJ is defined by the imag- 
inary part of (e d + e m ). 

Let us estimate the largest possible L^J . Below we 
will show that L^J defines the large- wavevector cut-off of 
the SPP dispersion law in the lossless metal case. Thus, 
while setting up experimental conditions for excitation 
of short-wavelength SPP, we only need to compensate 
the imaginary part of e m to the extent that the small- 
est SPP wavelength achievable according to eqs.(6,7) is 
of the order of L e JJ . After straightforward calculations 
this condition gives the largest effective 2D Planck length 
achievable in the experiment as 
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J 2D — 



{2 1 / 2 L e / f A ) 1/2 « 40nm, (11) 

where Ao is the free space photon wavelength at the 
frequency of the surface plasmon resonance. This value 
is clearly of the same order of magnitude as the SPP 
wavelengths observed in the experiment, which means 
that the effective Planck scale is an important parameter 
in quantum nanooptics. The value of L e JJ defines the 
thickness of the boundary layer of the nonlinear dielec- 
tric which exhibits the quantum highway mirage effect 13 . 
In this surface boundary layer the average refractive in- 
dex substantially exceeds the value of refractive index in 
the bulk (because of the Kerr shift caused by zero-point 
fluctuations of the SPP field) . In addition, this boundary 
layer experiences very strong quantum fluctuations of the 
dielectric constant: Ae d ~ 1—. As a result, surface plas- 
mons, which propagate through this fluctuating dielectric 
layer experience strong dephasing (loss of phase coher- 
ence) effect. In fact, according to Ioffe-Regel criterium, 
SPP with a wavelength shorter than L e JJ cannot propa- 
gate at all. Such SPPs experience Anderson localization. 
Thus, the value of Lf^J defines the high momentum cut- 
off of the SPP dispersion law in a way that is similar to 
the momentum cut-off at 1 / L p in quantum gravity^. 

Let us show that in many other respects the role of 
L^jJ parameter in quantum nonlinear optics is very sim- 
ilar to the role of the Planck length in quantum gravity. 
First, let us show that it defines the quantum uncertainty 
of the dielectric constant in a way that is similar to the 
quantum gravitational uncertainty relation for the space- 
time metrics. Expressed in terms of connection T, the 
quantum gravitational uncertainty relation is ArZ 3 > L 2 , 
while in terms of the metric tensor it may be expressed 
as 



*s=*"' e A,,,, - m A9l2 - L '- (12) 

where I is the linear size of a space-time regionitV 
As a result, the effective 2D Planck scale as seen by This uncertainty in the gravitational field comes from 
the surface plasmons the quantum fluctuations of space-time. 
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In the language of the effective metric (1), which is ex- 
perienced by SPPs in their propagation along the metal- 
dielectric interface, this uncertainty relation is translated 
into 



2Ae 



> 



L e Jf 2 
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(13) 



where Ltf* is obtained from eq.(5). It is easy to see 
that this expression coincides with eq.(8) from ref^, 
which was obtained via calculation of the zero-point 
fluctuations of the SPP field. Thus, the effective met- 
ric fluctuations experienced by SPPs are very strong: 
Agoo ~ goo, an d the dephasing effect of the dielectric 
constant fluctuations described above is a direct analogue 
of the phase coherence loss in quantum gravity^. 

Below we will show that similar to the Bekenstein- 
Hawking result for the black hole entropy^, the value 
of L^J defines the optical entropy bounds for a metal 
nanoparticle immersed in a nonlinear optical medium, 
and for a nonlinear dielectric microdroplet on a metal 
surface. As has been shown in refii^, the effective met- 
rics experienced by the SPPs near such objects may to 
some degree emulate the metrics of a real black hole. We 
will show that the entropy bounds in nonlinear quan- 
tum nanooptics may be expressed as the ratios of the 
droplet perimeter (nanoparticle area) to the effective 
Planck length (effective Planck length squared). These 
results reproduce the well known quantum gravity ex- 
pressions for the entropy of real black holes. 

Let us briefly recall the surface plasmon black hole 
analogy, which is described in detail in 18 . If a dielec- 
tric object (say, a droplet) is placed on the metal sur- 
face (see Fig.l), the SPP dispersion law will be a func- 
tion of the local thickness of the droplet and the local 
value of its dielectric constant, which may or may not 
be constant throughout the droplet. We may assume 
that the spatial distribution of the dielectric constant 
inside the droplet is chosen by an experimentalist such 
that at each illumination frequency in the range between 
oj p /(1 + erf) 1 / 2 and Wp/2 1 / 2 (these are the frequencies 
of the SP resonance at the metal-dielectric and metal- 
vacuum interfaces, respectively) there will be a closed 
linear boundary inside the droplet for which the SP reso- 
nance conditions are satisfied, and both phase and group 
velocity of surface plasmons tend to zero at this bound- 
ary. Quantitatively the effective metric experienced by 
surface plasmons near the droplet boundary may be writ- 
ten as 



z* 2 dt 2 



dx — dy 



(14) 



where c*(x,y) is the local SPP phase velocitjii^. The 
behavior of c*(x, y) is defined by the shape and thickness 
of the droplet near its edge, by the thickness of the metal 
film, and by the frequency of SPPs. In order to emulate 
the Rindler metric, the droplet geometry may be chosen 



so that c* — axe in the vicinity of x — 0. Naturally, 
only some features of real black holes are emulated in 
this model. The limitations of this model are described 
in detail ini&. 

This toy-model is also limited in terms of emulating 
only the motion of individual electromagnetic quanta in 
some static effective metric, which is unaffected by the 
motion and changing spatial distribution of other elec- 
tromagnetic quanta in the system. This limitation may 
be lifted if we consider an experiment performed with 
a droplet made of nonlinear dielectric material. If the 
optical nonlincarity of the droplet may be described by 
equation (2) with \^ > 0: the self-focusing of surface 
plasmons may lead to an effective "gravitational col- 
lapse" of the SPP field near the critical surface inside 
the droplet. As described above, this type of nonlin- 
earity causes an effective attractive interaction of SPPs 
with each other. Thus, we may imagine a situation in 
which a liquid droplet is illuminated with an intense plas- 
mon beam at a frequency below u> p / (1 + e^) 1 / 2 , so that 
a low intensity SPP field would not experience a criti- 
cal surface near the droplet edge. However, the increase 
in the droplet refractive index due to the high intensity 
SPP field will cause the appearance of the critical sur- 
face at which the surface plasmon resonance conditions 
are satisfied. As a result, the plasmon field will collapse 
towards this arising critical surface. In order to com- 
plete the picture, we should take into account the strong 
quantum metric fluctuations near the emerging critical 
surface. Similar to the case of metal-dielectric boundary 
described ini£, the average refractive index of this thin 
fluctuating layer (with the thickness of the order of L%jj) 
will be larger than the refractive index of the dielectric 
far from the critical surface. Thus, it would act as a self- 
induced waveguide and support SPP whispering gallery 
modes, which were experimentally observed iniS. 

Let us compute the optical part of the entropy S of 
such a toy SPP black hole/dielectric droplet produced as 
a result of the described SPP self-focusing experiment. 
We will do calculations in a usual way as S = InT, where 
r is the statistical sumi^. The dominant contribution 
to T is produced by SPPs near the frequency of surface 
plasmon resonance. For such SPPs the dispersion law 
is almost flat. While they differ in their wave vectors, 
all the SPPs have the same frequency. That is why a 
very large entropy state is created in experiments near 
u> sp . If we assume that we have N plasmons in such 
system, which are somehow distributed over M possible 
single-plasmon states in the boundary layer adjacent to 
the critical surface (see Fig.l), the statistical sum is given 
by 

_ (M + N - 1)! ^ M + N 1/a (M + N) M + N 

N\{M-l)\ ~ [ 2ttNM } N n M m 1 ' 

Since the experiment is conducted near the frequency 
of the SP resonance, the sum (ed+£ m ) is small. According 
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to eq.(2) the maximum number of plasmons in the system 
is limited by the value of 



Et 



-(e m + e d ) 

47TX (3) 



(16) 



(note that the value of td near the critical boundary is 
shifted due to zero-point fluctuations of the SPP field). 
Thus, necessarily M » N near the critical surface, and 
the statistical sum from eq.(15) may be simplified as T w 
(M/N) N . As a result, we obtain a simple expression for 
the entropy as S — Nln(M/N). The entropy bound is 
defined by the largest number of plasmons N max , which 
may be fitted into the boundary layer near the critical 
surface, and thus by E max from eq.(16). As a result, 
may be estimated as 



-2 

max 
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(17) 



and (neglecting the factor ln(M/N) and other numeri- 
cal factors of the order of one) we obtain S max ~ r/L^J , 



where r is the radius of the critical surface inside the 
droplet. In a similar fashion, the optical entropy bound 
for a metal nanoparticle immersed inside a nonlinear 



medium may be found as S n 



r 2 /L 



eff2 
2D ' 



where r is the 



radius of the nanoparticle. These estimates reproduce 
the well-known quantum gravitational limit on the en- 
tropy of a given spatial region established by Bekenstein 
and Hawking. While important for the field of quantum 
nonlinear nanooptics, these bounds represent an interest- 
ing example of a physical situation in which the entropy 
bounds have very transparent physical origins, and arc 
easy to calculate. 

In conclusion, optical entropy bounds of metal 
nanoparticles immersed in nonlinear optical media and 
nonlinear dielectric microdroplcts on metal surfaces have 
been calculated near the frequency of the surface plasmon 
resonance. Similar to the Bekenstein-Hawking result for 
the black hole entropy, the entropy bounds in nonlinear 
quantum nanooptics may be expressed as the ratios of 
the droplet perimeter (nanoparticle area) to the effective 
Planck length (effective Planck length squared). 
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